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Abstract 

In recent years, the paradigm of 'soft glassy matter' has been used to de- 
scribe diverse nonergodic materials exhibiting strong local disorder and slow 
mesoscopic rearrangement. As so far formulated, however, the resulting 'soft 
glassy rheology' (SGR) model treats the shear stress in isolation, effectively 
'scalarizing' the stress and strain rate tensors. Here we offer generalizations 
of the SGR model that combine its nontrivial aging and yield properties 
with a tensorial structure that can be specifically adapted, for example, to 
the description of fluid film assemblies or disordered foams. 

1 Introduction 

Many soft materials, including foams, dense emulsions, slurries, pastes, and tex- 
tured morphologies of liquid crystals, are characterized by the presence of struc- 
tural disorder on a mesoscopic scale (nanometres to microns), causing metasta- 
bility and slow dynamical evolution. Such materials are nonergodic, and can 
therefore be viewed, in at least one sense, as glasses. In the current work we 
address only shear-thin ning materials (for generalizations to shear thickening see 
Head et alJ fcoOll . 12002^ where data for the steady state flow curve are very often 



fitted either to (a) the 'Herschel-Bulkley' form, (7(7) — ay ~ 7 P , where ay > 
is a yield stress and < p < 1; or (b) to the 'power law fluid ' form which is the 
same except with ay = 0. (A survey of edible soft matter (jHolHsworthl . ll W± 



lists a very large number of instances of such fits in the literature.) The linear 
rheological spectra are also often close to power law fluid form, or else exhibit 
near constant storage modulus G'(uj) and anomalously flat or even rising loss 
modu lus G"(lo) as frequency is decreased. In many instances, aging effe c ts are 



seen 



dHohler et ail. Il99fll: ICloitre et all 1200(1 ICohen-Addad and Hohlerl . l200ll : 



Viasnoff and LecmeuxL 120021 : IViasnoff et al l . 1200.4 ICloitre et all 1200.4 ) in which 



the material gradually gets more elastic and less lossy as time goes by. These 
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effects can be quite complicated, with the rate of aging depending on stress 
( Cloitre et aT1 . l200(t l. 



2 The SGR Model 



Much of the aforementioned phenomenology, including some but not all aspects 
of the observed aging behavior, are captured by a simple and gener i c model 



called the SGR (or 'soft glassy rheology') model (|Sollich et all Il997l : ISollichl . 



1998: iFielding et all l2000t) . This model is base d on Bouchaud's trap model of 
glasses ( BouchaudT 19921 : Bouchaud et all Il995l ) , and envisages mesoscopic ele- 
ments whose dynamics consists of independent hopping between local traps (or 
free energy minima). In the context of e.g. a foam, such 'hopping' events cor- 
respond to yielding, where a cluster of bubbles rearranges into a new and more 
energetically favourable topological structure. The hopping is controlled by an 
effective temperature parameter x = T e g/T g which lies near a glass transition 
(x ~ 1). Such a transition exists if, as we shall assume, the distribution of trap 
depths is exponential. This choice allows many-body effects (which are undoubt- 
edly present near any glass transition) to be ignored without losing the transition 
altogether, and is therefore an intrinsic part of the trap model picture as usually 
formulated. 

In its original form (So llich et a l.. 1997; ISollichl . [l998) which addresses only 
simple shear flows, the SGR model combines the hopping dynamics of Bouchaud's 
model with the buildup of local elastic shear strains in the mesoscopic elements; 
these are assumed, between hops, to be linearly elastic and to deform affinely 
with the applied flow. Upon hopping, the local strain I is reset to zero. The 
stored elastic energy density kl 2 /2 (with k an elastic modulus) in each element is 
offset against the trap depth and lowers the local activation barrier to hopping; 
this leads to shear thinning. The dynamics of the SGR model is contained in the 
time evolution equation for the probability distribution P(E, I, t) for a mesoscopic 
element being in a trap of depth E (> 0) with local shear strain I: 



P(E,l,t) = -jdP/dl-T e 



-{E-H 2 /2)/x 



P(E,l,t)+T(t)p(E)5(l) (1) 



Here Tq is an intrinsic jump rate, and T(t) = J Toe" 



-(E-kl 2 /2)/x 



P(E,l,t)dEdl is 



the overall jump rate allowing for the modulation of barriers by the strain; p(E) = 
exp(—E) denotes the distribution of trap depths into which elements can jump. 
Note that in proper dimensional units, we should write (E - vkl 2 /2)/(xk B T g ) 
instead of (E — kl 2 /2)/x, with v the volume of an element. In the following, we 
choose to measure energy densities such as kl 2 /2 in units of k B T g /v, and energies 
such as E in units of k B T g , so that the extra dimensional factors disappear. The 
macroscopic shear stress is taken to obey 



a(t) = J P(E,l,t)kldEdl. 



(2) 



For a detailed discussion of how the a forementioned flow and aging phe - 
nomenology ar i ses fr om this SGR model, see ISollich et all dl997h : ISollichl |l99? 



Fielding et al. (|200f]l ). A review which puts this in a broader context is iGate; 
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(2003). There are, of course, many open issues with the model. One of these con- 
cerns the interpretation of the noise temperature x and whether or not this should 
depend on flow history: in this paper we assume it does no t. Also, the Theologi- 



cal aging predictions of the model, though surpris ing ly rich dFielding et allbond : 



Viasnoff and LeaueuxL 120021: IViasnoff et all p003l. do not include all those found 
experi mentally; see e.g. iHohler et al. (1999); Cloi tre et al. I (|2000r i:l Cohen- Addad and Hohlerl 



(2001). Nonetheless, the model represents a useful step towards understanding 



the r heology of materials that are not time-translation invariant (Fielding et al 
l2000h . 

A separate shortcoming of the SGR model lies in its tensorial simplicity. Be- 
cause simple shear flow and linear local elasticity are both assumed, no normal 
stresses can ever arise. The strain variable I in Eq^is effectively a scalar, as is 
the macroscopic shear stress a. The form chosen for the shear thinning is also 
quite restricted. In what follows we address these shortcomings by showing how 
the SGR model can be 'tensorialized' with minimal damage to the appealing phe- 
nomenology that it contains. This allows for various forms of nonlinear elasticity 
at the mesoscale, and also lets us consider arbitrary deformation histories rather 
than just simple shear. 

Our starting point is the constitut ive equation for the standard, scalar SGR 
model, whose derivation (jSollichl . 119981 ) we briefly recall here. We imagine that at 
the zero of time the sample is prepared in a known state that has, for simplicity, 
all mesoscopic elements unstrained; this assumption could be relaxed, but saves 
algebra. This state is characterized by P(E,l,0) = Pq(E)5(1). The constitutive 
equation then reads 



a(t) = G (s«,)fc7(t) + / nt')Gi{z tt ,) fc[ 7 (t) - 7 (t')] dt' 
Jo 



(3) 



Here Go(^o) is the fraction of elements present originally (at t = 0) surviving 
to time t > 0. The factor A; 7 (t) is the stress they contribute. The integral is 
over elements created at t 1 < t, with creation rate T(t') and survival probability 
Gi(zu') to time t; each has been strained through ^(t) — j(t') and con tributes 
stres s accordingly. The following equations complete the prescription dSollichl . 
ll99Sh : 

1 = G (V)+ / nt')Gi{z u i)dt' (4) 
Jo 



Ztt> 



t> 



exp 



(k [ 7 (t") -7(0] 2 /2x) dt" 



(5) 



G (z) 



(e-sroexpf-E/*)^ . G 1 {z) = {e- ZU]C ^- a/x) ) p . (6) 

Equation^Jfixes, by normalization, the jump rate T(t') required in Eq|31 Equation 
Eldefines an 'effective time interval' z between times t and t'; this is defined in such 
a way as to absorb the factor by which the jump rate out of a trap is enhanced 
due to the presence of a strain. In other words, within the scalar SGR model, 
the nonlinear effect of shear can be viewed as 'making the clock tick faster' for 
any element that has accumulated nonzero value of the local strain I since its 
creation. Indeed, this is the only nonlinearity in the model, which is why EqEl 



D -zTo enp(-E/x)\ 
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defines the survival functions for each class of element exactly as one would in 
Bouchaud's model, except that z replaces t. 



3 Tensorial SGR Models 

We now turn to the main agenda, which is to relax the tensorially naive assump- 
tions of the SGR model as formulated thus far. Fortunately the structure and 
interpretation of the scalarized SGR constitutive equation, Eq|3J allow one very 
easily to construct various tensorial descriptions based on the same underlying 
physics, but with more specific and fully tensorial models for the local elastic 
and yield behaviour of the mesoscopic elements. Because the connection with 
Bouchaud's trap model remains intact, these models will retain the interesting 
aging behavior and spectral phenomenology outlined in the introduction. The 
same applies to flow curves although these will be somewhat modified by any 
additional nonlinearities now introduced at the mesoscale. 

We present this first for a schematic dumbell based SGR model. A somewhat 
simpler class of models, closer in spirit to the scalar SGR picture, is then sug- 
gested. Members of this class adapted to describe the case of foams and dense 
emulsions, are then studied in more detail. 



3.1 A Dumbell SGR Model 

We consider an ensemble of dumbell-like objects. Each dumbell is characterized 
by an end-to-end vector u = u M of length u = |u|; the stress contributed by such 
dumbells is deemed to be = nk(u^u u ) p where the average is over a distribu- 
tion function P(E, u, t); n is the density of dumbells and k an elastic constant. 
(The units of k are, in this section, different from those used previously; e.g. 
kv? /2 has the dimensions of an energy and is therefore assumed to be expressed 
in dimensionless multiples of k^T g .) Thus far, apart from the appearance of the 
E variable, this formulation resembles a first step towards the upper convected 
Maxwell model ( Gates! . 200.^1 ) . where the dumbell dynamics consists of two par- 



ticles advected by the flow connected by a spring. (This is often used as a model 
for sub-entangled polymers.) In the present context we are not restricted to 
thinking of our dumbells as polymers; they represent unspecified elastic objects, 
each of which is held in place by its neighbors but can make discrete stochastic 
rearrangements that relax its stress locally. 

SGR-like dynamics (quite unlike that of the Maxwell model) is now introduced 
by assuming that each dumbell follows the flow afflnely, except that from time 
to time it makes a jump to a completely new configuration, with a jump rate 
roexp[— (E — ku 2 /2)/x]. Here E is an energy barrier which is lowered by the 
stored elastic energy ku 2 /2. (One can imagine the dumbells as hooked into a 
network of neighbors, which deforms afflnely; but when a particular dumbell 
becomes too elongated, its connections to the network will be more likely to 
break.) After a jump, the dumbell is assigned a new yield energy E drawn 
from the usual prior distribution p(E) = e . It is also assigned at random a 
new value of u, drawn from the equilibrium distribution at noise temperature x, 

Peq(u) ~ exp( — ku-1l/2x). 
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The resulting equation of motion, closely analogous to Eq^ is 



P(E,u,t) = -(dP/du)-K-u-T exp[-(E-ku 2 /2)/x]P + T(t)p(E) Peq (u) (7) 

with K the rate-of-strain tensor and T(t) = ro(exp[— (E — ku 2 /2)/x])p the total 
jump rate. A new feature is the appearance of j? e q( u ) m the last term to replace 
5(1) in the scalar SGR model: this is the closest we can get to the assumption, 
made there, of 'zero strain' in new mesoscopic elements without actually setting 
u to zero (in which case every dumbell would collapse to a point and never be 
stretched under the affine flow). This choice will recover Boltzmann equilibrium 
for the u (at temperature x) in the absence of flow, which is not true of the 
scalar model: but it is a m oot point wh e ther t his should be recovered, since 
x is not a true temperature ( Sollich et al. . 1997 : Sollichl . 1998t ). An alternative 



choice, probably not much different in practice, would be to use 47rp eq (u) = 
5 ^(u-u) 1 / 2 — (3x/k) l / 2 ^J , corresponding to the selection of new dumbells with 

random directions on a sphere of radius (u 2 )lq 2 = 3x/k. 

Now assume that at time zero a system is prepared by some definite process 
(e.g. a quench) that gives a known initial distribution of dumbell barrier heights 
and vectors P(E, u, 0) = Po(E)po(u) (chosen factorable for simplicity). Then, 
following the same arguments as lead to EqOl the constitutive equation can be 
written (with an additional coefficient a discussed below) 



a tr Mi/ (t) = (G (z t0 (u))nk( 

+ rr(t')(Gi(z ti ,(u))nA : (E^.u) At (E^-u) 1 ,) Pcq ^ / . (8) 
Jo 

Here E tt > is the deformation tensor between times t' and t, while Go,G\ are 
precisely as defined in EqEJ the total jump rate T(t) is found from 

1 = (G (z t0 (u)) P0 + f r(t')(G 1 (^(u))) Pcq ^ / (9) 



and the 'effective time' variable z, which is now explicitly dependent on the end- 
to-end vector u with which an element was created, obeys 

z tt ,(u) = y exp (A:|E^.u| 2 /2x) dt" . (10) 

The above constitutive equation involves, as well as a tensorial dependence 
on strain of the elastic stress and of the yield energy, an extra layer of averaging 
over the Uj variables describing the end-to-end vectors with which dumbells were 
created after their most recent jump. This makes its quantitative analysis (which 
must mainly be done numerically from this point onwards) rather cumbersome 
and we do not pursue it in this paper. Nonetheless, it is clear by the construction 
of the model that its behavior in simple shear flows will differ relatively slightly 
from the scalar SGR model; we should expect Herschel-Bulkley (x < 1) and 
power-law-fluid (x > 1) flow curves with power law G*(to) spectra at x > 1 and 
aging behavior, very similar to that described by Fielding et al. ( 2000h . for x < 1. 



But, unlike the scalar version, the model is capable of nontrivial predictions for 
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normal stresses under shear and also for the rheological and aging behavior in 
elongational and mixed flows. 



Th e model just presented has some features in common with that of Michel et al 
(|20Qlh . If the dumbells are thought of as polymer strands, it could be used to 
represent physical gels with cross links having a broad (in the model, exponential) 
distribution of activation energies. As formulated so far, these energy barriers 
are lowered in height by the full stored energy in a network strand between links. 
This is surely an exaggeration, since in practice only a part of that energy can be 
converted into the work of breaking a link. However, to improve this correspon- 
dence one could choose a > 1 in Eq|HJ In this case the stored energy appearing 
in the lowering of energy barriers is only a fraction 1/ot of the elastic energy of 
the strand, whereas the full value is used to calculate the macroscopic state of 
stress. (Equivalently, the plateau modulus in the model is now larger by a factor 
of a.) 

3.2 A Better Class of Model 

The extra averaging required above, in passing from the scalar to the dumbell 
SGR model, is rather a nuisance. Arguably, though, it ought to be redundant. 
After all, the basic elements of the SGR picture were not intended to repre- 
sent individial polymers or particles (as the dumbell idea tacitly assumes) but 
mesoscopic elements, within which a degree of local averaging can already be 
presumed. One could thus hope for a simpler description in which such elements 
have a distribution of yield energies E, but otherwise are taken as isotropic bodies 
with a well-behaved, deterministic elastic response, albeit in most cases nonlin- 
ear. No averaging over Uj would be necessary if we replaced our dumbells with 
elastic spheres, for example. 

The general structure of the SGR-type constitutive equation that results from 
this assumption of local averaging within a meso-element is as follows: 

a(t) = G (zto)Q(Eto) + / t r(t / )G 1 (z«OQ(E^)dt / (11) 

J o 

z tt , = f 'exp[R(E tH ,)/x]dt" (12) 
Jt' 

with Gq,G\ and T(t) obeying Eas l4l6l Here Q and R are tensor and scalar 
functions of E that can be freely chosen. Suitable choices for foams and dense 
emulsions are suggested below, but models for other specific materials could em- 
ploy quite different forms for these. By construction, the pre-averaged tensorial 
models defined by Eas lllll2*l are, for bland choices of Q and R, again expected 
to behave rather like the scalar SGR model (with glass transition at x = 1, aging 
and power-law fluid regimes etc.), while at the same time describing nontrivial 
normal stress effects under shear flow, and offering tractable models for the be- 
havior of soft glassy materials in extensional flow. Shear thickening models , along 
the lines developed from the scalar SGR model bv lHead et aP (120011 . l2002h . could 



also be introduced with suitable choices of Q and R; in particular, if Q is not 
strongly strain-thinning, a choice of R that drops sharply at intermediate strains 
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and rises again for l arge ones should lead to shear-thickening (and perhaps static 
jamming J behavior (|Head et all l200ll . l2002h . 



3.3 Tensorial SGR Models for Foams and Dense Emulsions 

Turning now to foams and dense emulsions, what are suitable choices for Q and 
R? The first of these determines (to within a factor of the plateau modulus G 
which we set equal to unity below) the stress created in a mesoscopic element 
as a result of a deformation. In particular, the instantaneous response to a 
nonlinear step strain is wholly controlled by Q. Therefore it fixes, for example, 
the ratio <p = N2/N1 of first normal stress differences in such an instantaneous 
response. Another significant quantity, closely related to Q, is the ratio x = G/J-'o 
between the elastic modulus and the stored free energy density J-q, which resides 
in the interfacial energy of the fluid films in an unstrained state. On dimensional 
grounds, G and T are proportional; each scales as <r/£ where a is an interfacial 
tension and £ a characteristic length scale associated with a foam droplet. 

Turning to the scalar quantity R, this represents the lowering of rearrange- 
ment barriers due to stored strain energy. In principle the lowering of a barrier 
could have a somewhat complicated dependence on the deformation E of the 
given element, but we assume for simplicity that it depends only on the free en- 
ergy density .F(E) in a mesoscopic element under strain. Moreover, again for the 
sake of simplicity, we assume this dependence is linear: 

R(E) = (13) 

where A is a parameter that should be, with our choice of units, of order one. 
(This parameter obviates the need for the parameter a introduced earlier in Eq|HJ 
it plays the same role of determining what fraction of the stored elastic energy 
can be used as work to overcome a rearrangement barrier.) With R defined as in 
Ea ll3l the ratio x can De read off from an expansion for small shear strain 7 as 
R = A X7 2 /2 + --, 

We now introduce three possible models for the elasticity of fluid film assem- 
blies (foams and dense emulsions) that predict approximate forms for both the 
step-strain response tensor Q(E) and the stored free energy J-(E) as a function of 
deformation. All expressions are given only for the relevant case of incompressible 
flows, det(E) = 1. 

Model 1. This model is inspired by the work of iDoi and Oht3 (fl99l) who 



addressed the rheology of emulsion droplets under flow (albeit with no attempt to 
address the foam limit, where droplets are in close proximity). For simplicity, Doi 
and Ohta assumed affine deformation of an isotropic assembly of fluid interfaces. 
Their result for Q may be written 

^'^T^/ 'pj^ ' (14) 

where the integral is over the surface of a unit sphere; the corresponding form for 
R (using EqllSjl is 

*<E) _ 1 f <f« , 



A Air J |E T .ul 4 
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Model 1 has \ = 4/15. For small shear strains one finds Q xx = (8/21)7 2 , Q yy = 
— (13/21)7 2 , Q zz = (5/21)7 2 and hence (p = —7/6, while for large strains (p = — 1 
as follows from Q xx = —Q yy /2 = Q zz = (5 /8)7. For unia x ial ex tension both Q 
and R can be found explicitly as shown bv lDoi and Ohtal dl99lh : see Eqs^MM 
below. 

Model 2. This is a restricted simplification of Model 1, again following 
Doi and Ohtal (|l99ll ) who presented a tractable analytic approximant for Q in 



the case of shear deformations only. In cartesian coordinates (x,y,z) with shear 
along x and gradient along y, the resulting Q for shear strain 7 is: 



Q( 7 ) = (l+7 2 /3)- 



-1/2 



/7 2 /3 

7 -2 7 2 /3 J (16) 
\ 




and the expression for the effect of strain on stored energy, Ea ll31 gives: 



A 



(l +7 2 /3 )l/2_ L ( 17 ) 



Model 2 has <p = — 1 at all strains and % = 1/3. 

Model 3. This is based on the work of Larson! ( 19971 ) who, unlike Doi and 



Ohta, was explicitly addressing the case of dense foams. He modelled the fact 
that, in such a foam, an affine deformation does not preserve the 120° contact 
angles between fluid films; it should therefore be followed by a very fast local 
relaxation which will not show up in rheology. Larson's model for the step strain 
tensor Q reads 

Q(E) = -2C 1 / 2 + ^(trC 1 / 2 )l (18) 

where C is the Cauchy strain tensor, obeying C _1 = E T E, and C 1 / 2 is the sym- 
metric tensor with C^-C 1 / 2 = C. The resulting stored energy expression gives 
(via EafTSl: 

R trC 1 / 2 

x = ^i— L (19) 

Mod el 3 has <p = — 3/4 at small strains and ip — > —1 at large ones; it predicts x = 
1/6 rtLarsonLll997ft. T hese a re considerably closer to the values found numerically 



( Eeinelt and KravniklEnnflh . at least for dry foams, than those for Models 1 and 2. 



This is hardly surprising since the Doi-Ohta analysis that underlies those models 
does not allow for the constancy of contact angles between films. (Nonetheless, 
for simplicity we stick to the latter f o r the numerics in the next section.) 

We note in passing that iLarsoul (|l997l ) suggests a number of other general 



forms for the tensor Q adapted, for example, to an assembly of elastic rather 
than fluid films; these could be chosen, within the tensorial SGR framework 
developed here, to address the flow of glassy assemblies of such objects. 

3.4 Numerical Results 

We now present a selection of numerical results, focussing on Model 1 and its 
close approximation for shear strain, Model 2. Throughout, we work in time 
units such that the microscopic time scale Tq 1 = 1. 
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The simplest case is a instantaneous step deformation E. We assume that 
before this is applied the system is in eq uilibrium, which impl ies x > 1 since for 
smaller x the SGR model exhibits aging rtFielding et all 1200(1) . From EqfTTl one 
then deduces, by arguments exactly analogous to those of Sollich (|1998l ). that the 
decay of the stress with time t after the step deformation is given by 



a(t) = Q(E)G eq (texp[R(E)/x 



(20) 



Here the time-dependence is given by the same equilibrium stress relaxation func- 
tion as in the scalar SGR model, 



G eq (z) 



(21) 



As in the scalar case, nonlinearities in the time-dependence arise only from the 
factor exp[R(E)/x], which produces a speed-up of the stress relaxation for large 
deformations. However, additional nonlinear effects occur in the tensorial model, 
through the instantaneous stress response Q(E). An interesting feature of Eo l20l 
is that all components of the stress tensor relax with the same time-dependence. 
This implies for example that the normal stress difference ratio </? = N2/N1 is 
independent of time, a prediction which could be tested experimentally. 

In figure n we show example results for step shear, obtained from Model 2. 
These demonstrate the nonlinearities in both the instantaneous stress response 
and the time dependence of the subsequent stress relaxation. Only the results for 
a xy are shown; the relaxation curves for N\ = —N2 only differ through a constant 
factor which reflects their different initial values. 

To emphasize that the tensorial SGR model can deal with deformations other 
than simple shear, we show in figure El analogous results for uniaxial extension, 
with deformation tensor 



(22) 



In this case we can use directly the unapproximated Model 1 since the in tegrals 
over u in Eas ll4ll31 can be performed analytically ( Doi and Qhtal . 1991 ). The 
stress tensor Q(E) is diagonal, with Q xx = Q yy = —2Q ZZ and 1 




Qz 



3 3e 



1) 



„2e 



+ 2e~ e + {e be - Ae 2t ) i{e 6t - 1) , r(a) = 



arctan sfa 



while the effect of deformation on stored energy is given by 

1 



R 



( - + e 2e i{e ie - 1) 

1 In lDoi and Qhtal Jl99lf) . the prefactor (e 3E — appears to have been omitted from the 
second term in the square brackets of Eg 1231 presumably due to a typographical error. For 
e < and hence a < 0, the function r(-) is defined by its natural analytic continuation, t(o) = 
artanh(^/— a) / yj—a. 



(23) 
(24) 
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Figure 1: Stress relaxation after step shear strain, calculated within Model 2 
at effective temperature x = 1.5. Strain amplitudes are 7 = 0.01 (circles), 0.1 
(squares), 1 (diamonds), 10 (triangles). The two plots show the effect of variation 
of the parameter A: from Eqs ll3l and!20 ( larger A implies greater speed-up of the 
stress relaxation. For 7 = 10 and A = 10 (right), this speed-up is so large that 
the stress relaxes to zero essentially instantaneously on the scale of the plot. The 
nonlinear dependence of the initial value of the stress on 7, which arises from the 
factor Q(E) in Eqj201 is independent of A and therefore the same in both plots. 
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Figure 2: Analogue of figure ^for step uniaxial extension. Shown is the relax- 
ation of a zz , calculated within Model 2 at effective temperature x = 1.5. Strain 
amplitudes are e = 0.01 (circles), 0.1 (squares), 1 (diamonds), 2 (triangles). For 
the smaller e, the initial stress response is linear, a zz = 2e, but at e = 1 and 2 
the expected nonlinear deviations from this become apparent. 
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Next we consider steady shear flow, with shear rate 7. From Ea llll one finds 
that the stress tensor in the steady state is given by 



/ 00 d7Q(7)G' 1 (7- 1 /o 7 ^ / e^ , )/-) 



(25) 



where Q(7) and are Q and R evaluated for a shear deformation with strain 7 
and 7', respectively. As in the scalar SGR model, one deduces that for sufficiently 
large x the integrals over 7 are dominated by values 7 oc 7 /To = 7 (recall that our 
underlying relaxation rate has been set to Tq = 1). This then gives Newtonian 
behaviour at small shear rates, with a xy oc 7 and N\, N2 oc -y 2 . As x approaches 
the glass transition at x = 1, however, values of 7 ^ 7 /To remain important for 
small 7. These cause non-Newtonian singularities in the low-7 behaviour, 



These power laws, and the yield stress behaviour for x < 1, are largely indepen- 
dent of the particular forms of Q and R. They only rely on a xy oc 7 and N\, 
N2 oc 7 2 for small 7, and on the growth of ^2(7) with 7 eventually limiting the 
largest element strains that occur in the steady state. An interesting observation 
is that the non-Newtonian effects in the normal stress differences manifest them- 
selves in a larger region above the glass transition (up to x = 3) than for the shear 
stress (up to x = 2). Another notable feature is that, when both are anomalous 
(x < 2), normal and shear stresses obey the same power law, in contrast to the 
familiar analytic case (a ~ 7, N ~ A f 2 ). Figure |3] shows some example results, 
calculated from Model 2 with A = 1. The crossover from yield stress behaviour 
for x < 1 to non-Newtonian power laws for x > 1 to Newtonian flow for x > 2 
or x > 3 can clearly be seen. For larger A, strain-induced yielding is more pro- 
nounced. This has two effects: it limits the element strains 7 in steady state, 
so that at given 7 both a xy and N\ decrease with increasing A. It also means 
that nonlinearities caused by strain-induced yielding are stronger, and deviations 
from the small-7 behaviour of Eg 1261 therefore appear for smaller 7. Numerical 
results for A = 10 (not shown) confirm these expectations. 

Finally, we consider a somewhat more complicated scenario, namely stress 
relaxation after cessation of a steady flow. From Eq llll one finds for this case 
that the stress at a time t after cessation of the flow is 



For t = 0, this just gives the steady state shear of EqEHl Ecis l25l27l generalize 
straightforwardly to other steady flows (and their cessation); e.g. for an exten- 
sional flow one merely replaces 7 and 7 by e and e everywhere. 

We now give a brief scaling analysis of the behaviour predicted by Eql271 
focussing on the regime where all timescales (7 -1 and t) are large compared to the 
microscopic timescale Fg , which equals unity in our chosen units. The reasoning 




for 1 < x < 2 
for x < 1 




for 1 < x < 3 
for x < 1 



(26) 




(27) 
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Figure 3: Results for steady shear, showing shear stress a xy (left) and normal 
stress differences (right) as a function of shear rate 7, as calculated within Model 
2 with A = 1. The curves from top to bottom are for a series of increasing noise 
temperatures x = 0.5, 0.6, ... 3, with x = 1 and x = 2 highlighted in bold. Note 
that, for small 7, the shear stress shows Newtonian behaviour <j xy oc 7 for x > 2, 
while for the normal stress differences non-Newtonian behaviour persists up to 
x = 3. For x < 1 the curves tend to nonzero limits for 7^0, demonstrating 
that the system exhibits a (dynamic) yield stress. 
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is similar to that used for other rheo logical predictions of the (scalar) SGR model 
( Sollichl . fl99^ : iFielding et all l200(t l . One uses the following facts. The function 



G\ has the initial value Gi(0) = 1; for z 3> 1 it decays as G\{z) ~ t~ x . The 
integral Jq dj' e R<a >l x is ~ 7 for small 7 but eventually grows very quickly with 
7, in fact typically exponentially since R^') grows at least as a power of 7'. In 
the denominator of Eg 1271 one can thus limit the 7- integration to the O(l) range 
before this exponential cutoff sets in. Separating the regimes where 7 < 0(7) 
(and thus the argument of G\ is 0(1)) and 7 > 0(7) (where the asymptotic 
power law decay of G\ is a good approximation), one finds that the denominator 
scales as 7 for x > 1 and as j x for x < 1. In the numerator of Eg 1271 one can use a 
similar analysis, though the small-7 regime does not need to be treated separately 
since for f > 1 the argument of G\ is always large. For -yt -C 1, the integral can 
again be cut off at 7 = O(l); depending on the value of x, it can be dominated by 
the regime of small 7 = 0(7) and one then needs to bear in mind that generically 
Qxy{l) ~ 7 while normal stress differences such as Q X x{"f) — Qyyil) scale as 7 2 
for small 7. For jt 3> 1, on the other hand, the cutoff in 7 is located where 
jt ~ Jq d"f' e R ^ ^ x and so typically grows logarithmically with jt; in the latter 
case strains 7 of 0(1) - up to the cutoff - always dominate the integral. Putting 
these elements together, one finds for the relaxation of the shear stress 



J 1 for t < 7 _1 

<Txy(t) _ 1 {jty x for t > 7 

ff*v(0) ~ J t 2 - x for t < 7 

\ t 2 ' x (^t)~ 2 for t > 7 

and for a typical normal stress difference such as N± 




(28) 



1 for t < 7 _1 

JVj(t) _ \ (jt)~ x for t > 7- 1 

iVi(O) ~ J t 3 ~ x for i < 7- 1 

t 3 -^^)- 3 for f > 7- 1 



for x < 3 
for x > 3 



(29) 



In Eas l28l2'9*l the power laws for t 3> 7 -1 are all subject to logarithmic corrections 
in 7t, which arise from the variation of the 7-cutoff discussed above. The scalings 
given apply to the generic tensorial SGR model of Eo llll including in particular 
our Models 1-3. The glassy nature of the model is manifested in the relaxations 
scaling with jt for low enough x, rather than with t as one would expect if the 
microscopic timescale dominates the dynamics. As in the results for steady shear, 
Eg 1261 glassy effects remain important up to higher x for normal stress differences 
(x = 3) than for the shear stress (x = 2). 

In figure |1] we show some numerical results, calculated from Eo l27l within 
Model 1 with A = 1. The power laws predicted by Eqs l28|2"9*l are quite well obeyed, 
subject to the expected logarithmic corrections for t 3> 7 -1 . As a general trend, 
we note that the normal stress difference N\ decays more slowly than the shear 
stress a xy . This is certainly true for x > 3, where iVj decays with a slower power 
law than a xy from Eqs l28l29l and for 2 < x < 3 where Ni remains essentially 
constant for t <C 7 -1 . From figure |1] we also observe the same trend for x < 2, 
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Figure 4: Results for cessation of steady shear flow, calculated within Model 2 
with A = 1, for four different noise temperatures x. Shear stress a xy and first 
normal stress difference N\ = —N2 are plotted as solid and dashed lines, respec- 
tively; in each graph a set of four curves are shown for each quantity, for shear 
rates 7 = 10 -4 , 1CP 3 , 1CT 2 , 1CT 1 from bottom to top. Theoretically predicted 
power laws from Eas 12812^1 are shown by bold dotted lines and annotated with 
their (negative) slopes. Top left: For x = 0.5; both a xy and N\ scale with jt. 
Logarithmic corrections to the predicted decay ~ (jt) ~ x are non-negligible in the 
range of jt shown. The underlying power law is nevertheless correct, as a fit of 
<J xy (t) for 7 = 10 -4 to the form [aln(7i) + b^jt)^ 1 / 2 shows (dashed-dotted line, 
just distinguishable around jt = 1). Bottom left: x = 1.5; the larger decay expo- 
nent (jt)~ x now makes logarithmic corrections less important. In the inset the 
curves are rescaled vertically by their value at t = to show that their time de- 
pendence is otherwise identical. Top right: x = 2.5. The shear stress a xy follows 
the predicted power laws from Eg 1281 quite well. In constrast to a xy , the normal 
stress difference only decays significantly for t 3> j , scaling throughout with 
jt. This is shown in the inset where (on the same vertical scale) N\ is plotted 
against 7^- Bottom right: for x = 3.5, both <J xy and Ni decay for t <C 7 
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however: even though Eas l28l29"l predict that a xy and JVi should have the same 
asymptotic decay shapes in this regime, N\ is seen to decay more slowly than 
a xy , presumably due to stronger logarithmic corrections. 

4 Conclusions 

In this work we have presented a general strategy for alleviating one of the more 
serious restrictions of the SGR model for soft glasses, namely its effectively scalar 
treatment of stress and strain, which results in an inability to consider normal 
stresses or flows other than simple shear. As might be expected, the tensorial- 
ization allows significant scope for tailoring the model to different classes of soft 
material, without altering the rather simple basic assumptions that underlie it. 
In particular, we have presented results for several tensorial SGR variants that 
offer a possible description of flowing foams and/or dense emulsions. Because of 
the amorphous packings and slow dynamics observed in these materials, they are 
prime candidates for the physical picture that underlies the SGR model. It is 
useful, we believe, to have a constitutive equation that can capture to reasonable 
accuracy the step-strain response of real foams (Model 3 is probably best in this 
respect) while also predicting nontrivial yield behavior, flow curves, and aging 
phenomena. These three features all arise in the glass phase of the model (x < 1) 
which is therefore the one most likely to be of interest in foam flows, except per- 
haps very close to the onset of rigidity which occu rs at volume fractions of the 
dispersed phase of around 58 % dMason et al.J.ll99fih 



As emphasized elsewhere dSollicli et all Il997l : Sollichl . 19981 : Fielding et al 



200GJ), the absence of tensorial structure is by no means the only shortcoming 
of the scalar SGR model; its various assumptions are all questionable at several 
levels and, even after tensorialization, the approach should not be viewed as com- 
plete in any sense. Experimental falsifications of its predictions are welcome, since 
these will help direct theoretical work towards more complete models. (Verifica- 
tions are, of course, also welcome.) Our development in this paper of tensorial 
versions of SGR should certainly allow more stringent comparisons to be drawn 
between theory and experiments. 

Indeed, in rheological terms one could argue that a proper tensorial treatment 
of stress and strain is the absolute minimum required for any kind of serious 
predictive modelling to begin. In this sense, the models presented in this paper 
could represent a 'coming of age' for the SGR approach to the rheology and 
rheological aging of soft materials. Despite the shortcomings of SGR, we are not 
aware of any competing approaches that directly confront the nonergodic features 
of these systems within a rheological context. 
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